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Abstract. - We extend the analysis of the mean field schematic model recently introduced 
Q for the description of glass forming liquids to the case of a supercooled fluid subjected to a 
shear flow of rate 7. After quenching the system to a low temperature T, a slow glassy regime 
is observed before stationarity is achieved at the characteristic time Tg. Tg is of the order of the 
usual equilibration time without shear r° for weak shear, 7r° < 1. For larger shear, 7r° > 1, 
local rearrangement of dense regions is instead enhanced by the flow, and Tg ~ l/(r7). 



Structural rearrangement in supercooled liquids and glassy systems is severely suppressed 
due to configurational restrictions, requiring a cooperative dynamics of correlated regions that 
involves many degrees of freedom |^ . This complex behaviour often results in a slow kinetics 
characterized by diverging relaxation times at the temperature of structural arrest To |^ and 
strong non-equilibrium effects, such as aging. A glassy system above To is generally observed 
to be off-equilibrium either because a modification of the control parameters, such as the 
pressure or the temperature, has been exerted or because it is driven mechanically. In the 
latter situation it is possible to show that the aging of some systems can be triggered by the 
external forcing and that they look younger if a larger drift is applied . This is also witnessed 
by the modalities of the violation of the fluctuation dissipation theorem ^. In many cases, 
by injecting power there is the possibility to interrupt aging and stabilize the system into 
a power-dependent stationary state. This effect is observed, for instance, in gelling systems 
under an applied shear force |Q. 

In this letter we study the out of equilibrium evolution of a glassy system, such as a 
supercooled fluid, quenched above To in the presence of a shear flow with rate 7 The 
analysis of the dynamics is carried out in the framework of a model [Q recently introduced 
for the description of the glassy behavior close to the dynamical transition. The approach 
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is simple enough to be handled analytically, allowing explicit calculations. Our main result 
is the observation of two regimes: a weak shear regime, where the flow does not practically 
affects the glassy behavior, and a strong shear situation, where the inner relaxation process 
is enhanced, because density fluctuations are convected by the flow, and aging is interrupted 
after a time Tg ~ l/(Tj). For t > Tg a shear-induced time-translational invariant steady state 
is entered where the decay of the two-time correlation function is characterized by the presence 
of inflection points. 

When a macroscopic motion of the fluid is present, the mean field version of the model 
introduced in is generalized by the following constitutive Equation 



dpjr, t) 
dt 



V • t)V{r)] = D{t)V^p{r, t) + r;(f, t) 



(1) 



where p{r,t) is a coarse grained particle density, and D{t) ~ (Af(p)) is the average mobility 
of the particles. 77 is a gaussianly distributed random field, representing thermal noise, with 
expectations (r?(r,i)) = and {ri{r,t)r]{r' ,t')) = -2TD(t)V'^ [5{r - f')5{t - t'% where (...) is 
the ensemble average and T is the temperature of the bath. In Eq. (|l|) the second term on the 
l.h.s. is an advection contribute due to the presence of the flow, which was not considered 
in [0; is the velocity of the fluid: for the case of a plane shear flow we consider V{r) = ^ycx, 
ex being the unitary vector in the flow direction. Eq. describes a system where the shear 
rate is homogeneous throughout the sample and applies therefore only to materials which can 
support such a flow (these are sometimes referred to as soft glassy materials, as opposed to 
hard glassy materials where strain localization or fractures are observed po|). 

The present model aims to describe the main features of the out-of-equilibrium dynamics 
above Tq. It is schematic in spirit and, in order to be generic, leaves aside as much system 
specific details as possible. The basic assumptions is that a good deal of the complex behavior 
of glassy systems can be encoded into the conventional convection-diffusion equation (|^) 
by means of a suitably chosen particle mobility M{p). In a quickly vanishing function of 
the density, M{p) = exjp{v[p — 1]^^} was proposed on phenomenological grounds; here w is a 
(temperature dependent) parameter and the particle density has been rescaled so that p ~ 1 
is the point of dynamical arrest. This form of the mobility has been obtained in different 
approximations by several authors in apparently heterogeneous contexts as the free-volume 
theory of the glass transition Q or in "car parking" problems in one dimension pl] |. In the 
mean field version of the model, expressed by Eq. ([^), M{p) is replaced by its average value 
that can be computed as 



D{t) ^ [2TrS^{t)] 



-1/2 



M{p} 



-(p-p)V[2S^(t)] 



dp 



(2) 



where p — (p) and S'^{t) = ((p — p)^), since the density distribution is gaussian 

From Eq. by transforming into momentum space, the following formal solution for the 

two-time correlator C(A:', A:, At) = {p{k' ,ti)p{k,ti+At)) = C{k,ti,ti+At)5 IC{At) + k' 

is found, with 



C{k,ti,ti + At) = x[/C(At),ti]exp 



At 



K?{s)D{ti +At-s)ds 



(3) 



where /C(s) — k -\- jkxSCy, Cy being the unitary vector in the shear direction. Notice the 



presence of the delta function 5 



lC{At) + k' 



distortion induced by the flow ||T2|. x(fc, t) 



, as opposed to the usual 5{k -t- /c'), due to the 
C{k,—k,t,t) is the structure factor that, with a 
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high temperature disordered initial condition x(fc,0) = A, evolves according to 

X(fc, <) - (A - r)e~' /o '^'(-)Dit-s)ds ^ y 
From Eq. (^) the average density fluctuations can be computed through 

'S'^(i) — {'^''^)~'^ ^\k\^{^x{k,t)dk, where A is a phenomenological momentum cutoff. For suf- 
ficiently long times the integral of the exponential terms in Eq. (^) can be extended to the 
whole /c-plane (letting K = oo) because the support of x shrinks towards the origin (see Fig. 1), 
yielding 

S\t) = ^(A - T)i?o(t)'^ • |i?o(t) - 7' 

where q = (S(i/d)[A/(27r)]'^, Yjd is the surface of the d-dimensional unitary hypersphere and 
Rn{t) = 2 z^D{t—z)dz. Eqs. (||,H) are closed coupled equations which allow the computation 
of D{t) and hence the whole evolution of the model, through the correlator (^). 

In the case with no flow it has been shown [|] that when p is close to the critical value 
p = 1 a cage effect produces a transient pinning phenomenon characterized by the constancy 
of the main observables, namely C^XiD and S. Then, for longer times, less dense regions 
start evolving, but yet high density regions are almost frozen and evolve slower, producing a 
glassy behavior characterized by the existence of many time scales and aging. The evolution 
of the system in this regime is characterized by a decay of C slower than the usual diffusive 
behavior with an exponential damping of the correlations. Aging is eventually interrupted after 
a characteristic time t° and an equilibrium state is entered where time translational invariance 
is recovered. For dense systems Jl3[ , equilibration is induced by thermal fluctuations which 
shake the frozen regions of the system. 

When shear is applied to the fluid a similar situation occurs. In the following we describe 
the main results of the analytical solution. We refer to a deep quench (low T) with a high 
density (p ~ 1). Initially the cage effect is observed, as discussed in 0|; this is reflected in 
Fig. 2 by the constancy of D for small times. For longer times the evolution starts and the 
glassy regime is entered, as for 7 = 0. The behavior of the structure factor in this time domain 
is shown in Fig.l where the anisotropic character of the correlations is evident. We consider 
the situations relative to different ranges of the strain ^t. 

Small strain. For sufficiently small ^t, neglecting the terms proportional to T and 7 in 
Eq. (||) one can solve Eqs. (||,||) obtaining 

D{t)^t-^{\ntY-^, (6) 

where 5 = 6/d, as without flow. This is shown in Fig. 2 where it is seen that the curves for 
D{t) initially collapse onto the 7 = line. Inserting the expression for D into Eq. (^) one 
finds that the terms proportional to 7 are negligible and C falls isotropically as an enhanced 
power law. Time translational invariance is lacking and the system ages. Shear induced effects 
become relevant at the crossover time Tc when the terms proportional to 7 and T cannot be 
neglected in Eq. (^). From the analysis of Eq. (||) we obtain 

r, =min{r;,7-^} (7) 

and r°, the time at which aging is interrupted in the corresponding undriven system, is 
computed in [Q. In the weak shear regime, with r° << 7"^, the glassy evolution is ended 
at f ~ r°, due to the thermal fiuctuations, and the presence of the flow only affects the 



Ro{t) 



-R2{t) 



(5) 
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asymptotic dynamics that will be discussed below. For the choice of parameters of Figs. 1,2 

9 



one has r° ~ 10® (this is the time at which D approaches the asymptotic constant value) and 



so, from Eq. (^, we expect the curves of D to depart from the 7 = case for 7 < 10~®, as 
shown in Fig. 2. 

Large strain. For large strain and values of 7 in a strong shear regime with 7^^ << r°, 
the glassy behavior is changed for t > Tc by the motion of the fluid. By keeping only terms 
proportional to 7 in Eq. (^) we find 

i?~e-(*/*°)''', (8) 

where is a constant. In this glassy regime the convergence to the asymptotic stationary 
state is enhanced by the shear. One can check that the term qT is negligible on the r.h.s. of 
Eq. (|) up to Tg = (T-f)-^, when aging is interrupted. The dependence of Tg on 7 is shown in 
the inset of Fig. 2, showing that Tg — Tg for small 7 while it approaches a, Tg ^ 7^^ law for 
larger values of the shear. This completes the description of the pre-asymptotic glassy stage. 

We consider now the behavior of the system in the asymptotic domain (t > Tg). Since 
Rn{t) 00, one has S{t) ~ qT so that D{t) — D{oo) — const.. This can be seen in 
Fig. 2. Hence one time quantities do not depend on t and two time observables, such as the 
correlator (^), are functions of the difference At alone 

C{k, At) = Texp{-D{oo)[2Atk^ + -f{At)^k^ky + ■^{At)^kl]} (9) 

Time translational invariance is then recovered. For 7 = 0, from Eq. (^, C has a simple 
exponential form. When the shear is present, on the other hand, the decay is faster for large 
At (except at k^ =0), indicating that the relaxation of fluctuations is enhanced by the flow. 
Actually the behavior of C is more complex, as shown in Fig. 3. In this picture the behaviour 
of C{k,At) is plotted at 7 = 10"^ for different k. In the sectors k^ky > all the coefficients 
of the polynomial form in the argument of the exponential in Eq. (^ are positive. Then C is 
depressed with respect to the case 7 = for any At. On the other hand, C develops inflections 
in the sectors with k^ky < 0. From the comparison with the case 7 = one sees that this 
effect causes correlations to be initially (for -fAt < —3ky/kx) enhanced by the shear and only 
successively dumped. 

Integrating the expression of Eq. (^ over the space of wave vectors we obtain the behaviour 
of the autocorrelation function 

f dk ^ T 

A{At)^ ——C{k,At) = 5 (10) 

A.<A(27r)'i 2¥7rf [i?(oo)Ai]* [l + ^(At)V] 

For At < At* — 4/(V77) shear is ineffective and A decays as for 7 = 0. At large strain the 
term proportional to 7 in Eq. (^0|) prevails so that A ~ 'y~^{At)~^-^^'^^'^K Correlations are 
then suppressed by the flow for At > At* . A relaxation time proportional to 7"^ is also found 
in molecular dynamics simulations of supercooled liquids in shear flow fl^ . 

In this paper we have studied the behavior of a quenched glassy system under the effect of 
an applied shear flow. The analysis has been carried out analytically in the framework of the 
mean field model recently introduced in Two shear regimes can be distinguished which 
behave differently. While for weak shear the off-equilibrium evolution is unaffected by the flow, 
for strong shear the evolution towards the stationary state is enhanced and aging is interrupted 
after a time Tg = l/{Tj) by the flow itself. We have also described how the properties of the 
time correlation function in the stationary state are modified by the flow. These predictions 
are amenable of experimental checks; it would be interesting to know if similar feature can be 
observed in real glassy systems such as hard spheres or colloidal suspensions. 
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Fig. 1. - Four snapshots of the structure factor x{k, t) are shown for d = 2 at times t = 2,7, 20, 60 for 
p = 0.95, T = 10"'' and 7 = 1. kx increases from left to right (fcx = is in the middle); ky increases 
towards the upper part of the figure. 
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Fig. 2. - D{t) is shown in d = 2 for p = 0.95, T = 10"'' and different values of 7 (7 = 
1, 10"SlO"^,10"^,10"®,0) from left to right). The curves with 7 = 10"* and 7 = are not 
distinguishable in this plot. The finite final value D{oo) is very small for high densities and practically 
coincides with the t axis in this picture. In the inset the characteristic time Tg when aging is interrupted 
is plotted against 7. The solid line is the 7"^ law. 
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Fig. 3. - The two-time correlator C{k, At) in the stationary state is plotted against At for 7 = 10~^, 
ky = 1 and different values of kx and kz as indicated in the legend. Solid lines refer to wavevectors in 
the sector with kxky < whereas dashed lines correspond to the region kxky > 0. 



